Hard gluon damping in hot QCD 



A. Peshier* 

Institut fiir Theoretische Physik, Universitdt Giessen, 35392 Giessen, Germany 

The gluon collisional width in hot QCD plasmas is discussed with emphasis on 
temperatures near Tc, where the coupling is large. Considering its effect on the 
entropy, which is known from lattice calculations, it is argued that the width, which in 
the perturbative limit is given by 7 ~ 5^ ln{l/g) T, should be sizeable at intermediate 
I temperatures but has to be small close to Tc- Implications of these results for several 

I phenomenologically relevant quantities, such as the energy loss of hard jets, are 

pointed out. 



(N 



O 



I. INTRODUCTION 



^ ■ The dispersion relation and the damping rate of single-particle excitations in many- 

^ ■ particle systems are closely related to a variety of phenomenologically important quan- 
■ titles. In a QCD plasma at temperatures much higher than the transition temperature 
^ . Tc ~ 200 MeV, where the coupling g is small and perturbation theory should be applica- 
. ble, the quark and gluon excitation energies follow directly from the real part of the 1-loop 
^ ! self-energies, which are of the order {gTy. The calculation of their width, however, requires 
'-^ ! already to lowest order in g a summation of infinitely many diagrams. Resumming hard 
I I thermal loops (HTL), Braaten and Pisarski [1] obtained the widths of quarks and gluons 
D I at rest, which are proportional to g^T. The case of excitations with a finite momentum 
^ ' is more intricate because their width diverges even after the HTL resummation due to the 
exchange of soft magnetic gluons. With a cutoff of the order of g^T, either a magnetic mass 
^ ■ and/or the width itself, the width of moving charged excitations is ~ g^ln{l/g) T on rather 
^ ■ general grounds [2]. Even in the case of weak coupling, the gluon width tests directly the 
nonperturbative sector of QCD. 

A similar breakdown of perturbation theory occurs also in the calculation of the thermo- 
dynamic potential ^{T) = —p{T)V, at order 0{g^) [3]. The expansion in g, which is known 
up to 0{g^ lng~^) [4], is not reliable in the physically interesting regime as probed in present 
relativistic heavy ion experiments. For the large coupling strength expected at temperatures 
near Tc, it does not converge but shows a behavior typical of asymptotic series. In fact, one 
can hardly expect a converging expansion since it has to be defined in a circle in the complex 
plane, while in QCD a transition (7^ — —g"^ is presumably non-analytic. In perturbation 
theory, this manifests itself in the number of diagrams increasing rapidly with the order. A 
strategy to remedy the situation in practical calculations is a partial resummation of the per- 
turbative series, taking into account those classes of diagrams whose number increase most 
rapidly. These are related to the various orders in the loop expansion in the ^-derivable 
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approximation scheme [5] to be utilized below. In this scheme, the thermodynamic potential 
is expressed in terms of dressed propagators, which are determined sclf-consistently. The 
restriction to the leading loop order (s) for large coupling may seem counter-intuitive. Since 
it is essential for the following, it is worth mentioning another, yet related, motivation: An 
asymptotic series, to give the best approximation possible, should be truncated at an order 
related inversely to the coupling; for the QCD thermodynamic potential near Tc possibly 
already at the order 0{g'^) [6]. At large coupling, such a perturbative result is, however, not 
thermodynamically consistent since various thermodynamic quantities are related to each 
other by derivatives with respect to the temperature. Since T is also the relevant scale in 
the running coupling, which is large, a thermodynamically consistent approximation (when 
expanded in g) has to contain some contributions of higher orders. 

A truncation of a resummation scheme based on 2-point functions is, a priory, delicate for 
QCD because of gauge invariance. This problem can be evaded by receding to approximately 
self-consistent resummations of the thermodynamic potential using relevant gauge-invariant 
contributions to the propagators. Indeed, results calculated with HTL propagators [7, 8] 
agree with QCD lattice data down to temperatures of about 3Tc [31]. The HTL propagators 
can be reduced even further by neglecting the Landau-damping parts and retaining only the 
dominant pole contributions, approximating as well the dispersion relations by the asymp- 
totic mass shells. The resulting phenomenological models [10] can describe the lattice data 
even close to Tc because they allow for an IR-enhancement of the running coupling com- 
pared to the 2- loop formula used in the HTL calculations [32]. In all of these approaches 
the observed decrease in the effective degrees of freedom near Tc is directly related to the 
temperature dependent mass scale m ~ gT that characterizes the excitations. Interpreted 
as quasiparticlcs, they become heavy near T- due to the running coupling. While apparently 
this reflects important interaction effects (as motivated above), so far none of the approaches 
takes into consideration the expected width of the quasiparticles. This, however, is a priori 
not justified for large coupling, when the width might become comparable to the mass of 
the quasiparticles [12]. 

In principle, the dressed propagators and the widths can be calculated, by Schwinger- 
Dyson equations, in the ^-derivable approximation scheme. However, apart from the afore- 
mentioned sensitivity of the width to the soft QCD sector, there is the basic requirement 
of gauge invariance of physical quantities such as the width itself or the deduced pressure. 
Moreover, the resummed propagators need to be renormalized nonperturbatively. Notwith- 
standing the recent progress in these issues [13-15], the problem is involved. It therefore 
seems interesting to ask a reversed and simpler question: What can be inferred about the 
propagators from other quantities which can be reliably calculated by other means? An ev- 
ident choice for such a quantity is the thermodynamic potential, which can be numerically 
calculated rather precisely on the lattice. 

Clearly, one cannot expect detailed information from the bulk properties of the many- 
particle system, but estimates for some important features seem feasible. Indeed, the large 
quasiparticle masses 'predicted' in [7, 8, 10] compare well with direct results from lattice 
QCD [16]. It is obvious from the available phase space, though, that only properties of 
excitations with hard momenta, k ^ T, can be accessible. In the following I focus on the 
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entropy, s = —dft/dT. Since it provides a measure of the population of phase space, one 
expects an increased entropy for a system of off-shell particles as described by the width. 
From the estimate 7 ~ f?^ ln(l/5') T, the width might become large with increasing g, leading 
to the principal question whether this would be reconcilable with the small entropy near T^, 
as calculated on the lattice. 

To approach this question, in this paper some general relations between propagators 
and entropy will be discussed. Section II starts with a brief outline of the formalism of 
self-consistent approximations. In Section III, the case of particles with a Lorentz spectral 
function is considered in some detail, followed by an analysis of the sensitivity of the results 
on the spectral function. For the sake of transparency a scalar field theory is discussed before 
switching over to QCD in Section IV. Given the remarkably universal scaling behavior of the 
QCD entropy for various numbers of quark flavors [17], I consider here the representative 
case of the quenched limit of QCD. In the conclusions, some implications of the present 
findings are pointed out. Some formal details were deferred to the Appendix. 



II. PROPAGATOR AND THERMODYNAMICS 

Following the work of Luttinger and Ward [18], the thermodynamic potential of a 
system of particles with a given interaction can be expressed in terms of the exact 2-point 
function(s). Considering for simplicity a scalar theory with the propagator A, the expression 
reads in the imaginary-time formalism (setting the volume V — 1) 

(in(-A-i) + ha) - $[A] . (1) 

The self-energy 11 = Aq^ — A~^ is given by the sum of all 2-particle irreducible skeleton 
diagrams, which is closely related to the expansion of the functional $. The right hand side 
of (1) is a functional to be evaluated with the exact propagator A, which is obtained from 
the stationarity condition 5Q[A]/5A = 0. The functional variation leads to 

n = 2^, (2) 

i.e., the self-energy is self-consistently obtained by cutting a full propagator line in the 2- 
particle irreducible skeleton expansion of For a («;/3!)0^ -|- interaction, the 

functional reads 

$ = 30+3 QO + 12 Q + ■ ■ ■ 

From this exact representation, symmetry-conserving ('<I>-derivable') approximations [5] fol- 
low by truncating the ^-expansion at a given loop order. For a given order, the self-energy 
is then calculated self-consistently according to Eq. (2). 

To derive the entropy in terms of the resummed propagator, the Matsubara-sum in (1) 
is first transformed into a contour integral in the complex energy plane. After wrapping the 
contour around the real axis one obtains first 

Q = n{uj) Im(^ln(-A-^) + HA^ - $ , (3) 
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where A now denotes the retarded propagator; = Ji^a J duj/(2n), J^s = J d^k/ (2n)^, and 
n{Lj) = (exp(u;/T) — is the Bose- Einstein distribution function. Taking Sfl/SA — into 
account leads to 



on f On „.\ (9$ 

s — - 



(4) 



dT Jk^dT \ ' ' J dT 
This expression can be decomposed, s = + s', where 

.dip ^ _ ^ (i^ ln(- A-^) + Imn ReA^ , (5) 



s 



and 

, r dn ^ ^ 
s'^- — RenimA+ — 



(6) 

A 

Remarkably, in s' the first term cancels the contributions from $ with one and two vertices 
shown diagrammatically above, see e.g. [7]. In other words, the contribution (5), given the 
propagator, is the leading-loop resummed entropy in the ^-derivable scheme. According to 
the arguments put forward in the Introduction, it is a preferable approximation of the exact 
entropy at large coupling [33], 

s ^ s'^P . (7) 

The cancellation in s' is basically a topological feature, and expressions analogous to (5) 
hold also in other theories [7, 19], including QCD. In the context of the Fermi liquid theory 
it is called the dynamical quasiparticle contribution to the entropy [20]. Although the ap- 
proximation (5) is a nonperturbative resummation, it has a simple 1-loop structure and it 
does not depend on the vertices. The integrand can be rewritten using 

Imln(— A~^) — 7rsgn(lmA) — arg(A) 

= 7rsgn(ImA) e(ReA) - arctan(ImA/ReA) . 

In the first term, the real part of the propagator is negative for small a; > 0, as shown 
below, and it changes sign at o;^ = col- Using sgn(ImA(a;)) = — sgn(a;), this term yields the 
expression of the entropy of free bosons with the dispersion relation cuk (and zero width). 



i J^^ {-T In (l - e-^^/^) + CO, n{u;,)) . (8) 



Below, ujk will be simplify referred to as the dispersion relation of the dynamical quasiparti- 
cles although it is not the real part of the pole of the propagator (unless its imaginary part 
vanishes). In the total entropy, 

s = + As , (9) 

(where now the relation (7) simplifies the notation) the effects of a non-zero spectral width 
are solely due to the second contribution. 
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where A = ImA/ReA. For later reference it is noted that the second term in the parenthesis 
is ImA times the derivative of the phase arctanA with respect to ImA. 

The expectation that the entropy is increased for a non-zero spectral width, As > 0, can 
be verified from the representation (10) under rather general assumptions. To this end, the 
analytic propagator, for complex ko, is first expressed in the Lehmann representation, 

J-oo 2tT Kq — U! 

The spectral function is the discontinuity of the propagator at the real axis, 

p{u, k) = A''{u - ie, k) - A"(a; + ie, k) 
It is real, odd in cu with ujp{ijj) > 0, and it satisfies the sum rule 

^^P('^,fc) = 1 

for all values of k. This implies that the propagator approaches the free limit at large /cq, 

A«(..,fc)=r^^.|i^^ (13) 

Jo T: Kq — UJ^ Kq 

These general properties of the spectral function have several implications for the retarded 
propagator A(u;) = A"(a; + ie). Its imaginary part, which by the refiection principle is 
— I p(a;), satisfies 

ImA(a; = 0) = 0, 
ImA < for u; > 0, 
ImA — > for ou — > oo. 

Similarly, one readily infers 

ReA(cu = 0) 
ReA u}-^ 

and that odd-order derivatives of ReA(a;) vanish at a; = 0. Now consider a generic spectral 
function with a prominent peak near cuk and a characteristic width 7, possibly with some 
additional minor substructures. In this case ReA(a;) changes its sign only once for a; > 0, 
i. e., the 'dispersion relation' Uk is unique, which will be the only assumption for the following 
argument. In principle then, there are two typical cases of propagators, see Fig. 1: (i) the 
imaginary part is regular, and the real part is smooth, and (ii) the imaginary part is singular 
(but integrable due to (12)), and the real part is discontinuous. Common to both cases is 
that LOk is determined by the real part of the self-energy, 

ReA-\u;k) = AQ^Uk) - ReU{u;k) =0. 
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(12) 



/ 



du' p{uj') 



27r uj' 
for UJ — > 00, 



< 0, 
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FIG. 1: The real and the imaginary parts of the propagators and Api defined in Eq. (16) and 
in Sec. HIE, respectively. All quantities are in units of the width 7, and the energy scale E in the 
spectral functions is chosen such that cjfc = \/2. 




FIG. 2: The functions A = IniA/ReA, and /(A) = arctan A — A/(l + A^) (occurring in Eq. (10) in 
the integrand of As) for the propagators shown in Fig. 1. 

For singular spectral functions the dispersion relation cj^ and the peak position in p coincide, 
while they can be separated by an energy of the order of 7 for regular spectral functions. 
The integrand of the entropy contribution As is discontinuous at tu^; shown in Fig. 2 is the 
parenthesis term in Eq. (10). This factor is, in an approximate way, symmetric near tUfc, 
where A(ci;) = ImA/ReA is singular. Since the second factor in the integrand (10), 

dn _ u; exp{uj/T) 

dT (exp(cu/T) - 1)2 ' ^ ^ 



is monotonically decreasing with 00, it is plausible that As > 0. A rigorous argument is 
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given in the Appendix. From Fig. 2 it can also be expected that the entropy increase is 
smaller for singular spectral functions than for regular ones. 



III. SPECTRAL FUNCTIONS 

A. Lorentz spectral function 

In the previous section, the spectral function was introduced to deduce some general 
properties of the propagator and the entropy. At the same time, lacking a direct calculation, 
it is more intuitive (and more efficient due to the analytic properties) to model the spectral 
function rather than the propagator. From the spectral function of free relativistic particles 
with Aq ^ = /cq — o;^ where a;^ = + fc^. 



an often used ansatz to describe non-zero width is obtained by replacing the 5-function by 
a Lorentzian, 

^'^""^ = E [{u;-Ey + ^^ ~ J^JTW^) ■ ^^^^ 
The corresponding retarded propagator can be easily calculated by a contour integration. 



In general, the analytical continuation of the retarded propagator to complex energies is 
analytic in the upper plane. In the present case it has poles in the lower plane, at kg — 
±E — i'j. The parameter E is directly related to the dispersion relation. Choosing E'^{k) — 
— 7^, the propagator becomes 

A4^,fc)= ^ ^ . (16) 

With this convention [34], the parameter corresponds directly to the real part of the 
retarded self-energy. This has the advantage that the dispersion relation does not depend 

on 7, LOk = UJm- 

Turning now to the entropy, one should note that in general the mass and the width 
parameters are momentum-dependent. The resulting effects will be considered below; for 
now the parameters arc assumed to be constant. For the propagator (16), the contribution 
(8) to the dynamical quasiparticle entropy. 



,(0) 



(m) = ^ j^^ (-rin(l - e--/^) +a;^n(u;^/r)) , (17) 



is simply the entropy of free bosons with mass m. Corresponding expressions for QCD have 
been the starting point in the approaches [10], which interpreted the thermodynamically 
relevant transverse gluon and quark particle-excitations as quasiparticles with masses given 
by the asymptotic self-energies (and respective degeneracies) . 
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The contribution (10) due to the non-zero width reads exphcitly 



AsL(m,7) 



r dn I 
Jk^df [ 



arctan ■ 



27C0' 



— 2'~fuj- 



(18) 



^2 - a;^)2 + (270;) V ■ 

A numerical integration shows - in hne with the general expectation - that the total entropy 
Sl = sf} + Asl increases with the width and decreases with m, cf. Fig. 3. An interesting 
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FIG. 3: Left: The entropy siim, 7) as a function of 7 for several values of m (7 and m are in units 
of T). The dotted lines show the expansion (19) for m = and m/T = 1. Right: Contour plot of 
sl/sq', the contour spacing is 0.25, and the straight line marks s{m,^) = sq. 



detail is that SL{rn = ■y) is equal to the Stefan-Boltzmann entropy of the massless ideal gas, 
So = ^ TT^T'^. This is proven in the Appendix, where also the expansion 



SL(m,7) = So 



15 

7^ + 



15 72 



87r2 T2 87r2 



(19) 



for small values of m and 7 is derived. It is interesting to note that this result can be obtained 
by the expansion of the contribution s^^ with complex masses, SL(m, 7) ~ h{sL\m + i'y) + 



(m — 27)). 



B. Momentum-dependent mass and width parameters 

Due to phase space, thermodynamic bulk properties are mostly determined by hard mo- 
menta. Therefore, the entropy is expected to be not very sensitive on the exact momentum 
dependence of the width as well as on the dispersion relation (described by a momentum 
dependent mass parameter) at soft momenta. To quantify this expectation, the squared 
mass and the width are varied for < T by some factor from and 7, which are now 
considered as the asymptotic values. Denoting the resulting entropy by s^, the quantity 
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provides a measure of the momentum sensitivity. As shown in Fig. 4, r is indeed only of 
the order of a few percent when varying the dispersion relation. The sensitivity to the low- 
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FIG. 4: The sensitivity of the entropy sl on the momentum dependence of (left) and 7 (right). 
The parameters are varied for k < T (sec text) in the range [3,4] (outer band) and [^,2] (inner 
band). In both cases the asymptotic mass is m = 4T. 



momentum behavior of the width is even less. The Figures show the sensitivity for a rather 
large asymptotic mass; for smaller masses the sensitivity is even lower. This quantifies the 
expectation that the entropy is, to the extent required below, insensitive to details of the 
propagator at soft momenta. 



C. Specific shape of spectral function 

For a first test of the sensitivity to the specific form of the spectral function, let's consider 
the normalized function 

which has a more pronounced peak than the Lorentzian (15). It can be expressed in terms 
of the function with complex width parameters, 

Pq[i] = (VipL[Vij] + ^/^pL[V^J]) . (21) 

An analogous relation easily allows to obtain the corresponding propagator in terms of the 
Lorentz propagator (16). Replacing furthermore 7 — > \/27, the result reads 

^ , , , + 2ab'^ - 4ib^ 
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where a = u'^ — u'^ and b = 2'yu. The propagator (22) coincides with Ai{uj, k) for c<j ^ and 
uj — > ±oo as well as on the common mass shell = u"^. Thus, differences in the entropies 
can indeed be attributed to the spectral form rather than to a change in the dispersion 
relation. 

The differences in the propagators Al and Aq, apart from the analytic structure, are 
lucid for typical values of the 4-momentum, see Fig. 5. Nonetheless, there is almost no 
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FIG. 5: Left: The spectral functions of the propagators Al and Aq with u>m = (all quantities 
are in units of 7); if 7 = T, this corresponds, e.g., to k = m = T. Right: The corresponding 
entropies as functions of 7 and m (in units of T); for small masses the lines practically coincide. 

effect on the entropy even for large values of m and 7. This example leads to the question 
which features of p the entropy is actuahy sensitive to. Already expected in Sec. II were 
differences between the two generic types of spectral functions shown in Fig. 1. In any 
case, the integrand of the contribution (10) is notably different from zero only for energies 
— cijjtl ^7, cf. Fig. 2. A way to focus on this relevant interval is to consider large times, 
t > 7~^, in the Fourier transform p{t) of the spectral function. 



D. Spectral functions in Fourier space 



The Fourier transform of the spectral function is defined by 

duj 



Pit) 



-iujt 



27r 



duj 



Auit 



The sum rule (12), 

1= / —up{u;)= / dtp{t) / —ue'' 
which after a partial integration becomes i jdt f)(t)5{t)^ then translates into 

dp{t) 



dt 



t=o 



(23) 



(24) 
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It is plausible that the sum rule tests the small-i behavior in Fourier space since in momentum 
space it is closely related to the fact that the propagator approaches the free limit at large 
energies. The Lorentzian spectral function (15), with 

sin Et 
PL{t) = exp(-7|t|) — 

obviously comphes with the condition (24). Prom the fact that the oscillations and the 
attenuation of p{t) are related to the position and the width, respectively, of the peak of 
p(a;), one can easily construct other conceivable spectral functions. In the ansatz 

, , ^, , sin Et 
pAt) = f{t)—^, (25) 

functions with /(O) = 1 provide candidates for possible spectral functions if a;p(a;) > is 
satisfied. 

The similar entropies for the propagators and Aq, cf. Fig 5, may now be attributed 
to the similar large-t behavior of the spectral functions (basically also PQ{t) decreases expo- 
nentially, see Eq. (21)). Before studying this in more detail, it is noted that in terms of p{t) 
the retarded propagator reads 

roo 

A(A;o)=rM dt^^^'pit). (26) 

For the ansatz (25), 

= 2eL " ^^^^ 

= ^ {Hko + E)- - E)) , (27) 
the propagator can be expressed by a Fourier transform of the function /(i), 

J^icu) = / dte'^^'fit). (28) 



E. Non-exponential time behavior 



Although often assumed, an exponential decrease of p{t) is not dictated by any funda- 
mental requirement [21]. As already argued, the contribution (10) to the entropy, due to 
the non-zero width, is determined by the long-time behavior of the spectral function; it will 

increase if p{t) decreases faster, either by a larger value of 7 or due to the functional form of 
p{t). In the following, this is demonstrated systematically by some models for the spectral 
function (25), which are summarized in Table 1. 

In the models Pi and P^, pit) ~ fit) decreases asymptotically as t~^. This implies that 
the spectral function 



p{u) = -2ImA(a;) = {J'cioj + E) - Tj^u: - E)) , 



(29) 
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model 
L 
Pi 
P2 



Pn 

Pt 

G 



fit) 
exp(-7|t|) 

(l+7|t|)-« 
(l + (7t)2)-V2 

exp (-(7t)') 





e-*^r(0,-zx) 
1 + i2;e-*^r(0, -ix) 

(l+ix[7^„_i])/(n-l) 
iro(|x|) + fsgn(x)[/o(kl)-Lo(|x|)] 
i^^exp (-(x/2)2) [1 + Erf(ix/2)] 



TABLE I: Damping models (25) for together with their Fourier transforms J- defined in 

Eq. (28). r(0,2;) denotes the incomplete T-function, Kq is the modified Bessel function of the 
second kind, Lq is the modified Struve function, Erf is the error function, and x = lo/^. 



where J-'df^) = dt cos ut f (t) is the cosine transform of f(t), diverges logarithmically at 
uj = ±E, see Figs. 1 and 6. For 7 — 0, the free propagator is recovered. For the model 
Pi, e. g., this follows from r(0, x) = e^{x~^ + . . .) for x 00. Moreover, the incomplete F- 
function F(0,2;) is discontinuous at the negative axis. Accordingly, the retarded propagator 
A PI has cuts in the lower ko-plane, starting at the singularities. The entropies for the models 
Pi and are shown in Fig. 6. As anticipated before, the entropy increases with 7 much 





FIG. 6: Comparison of the models Pi and Pi with L anafogous to Fig. 5. The entropies almost 
coincide for the cases Pi and Pf. 

less than in the case of a regular spectral function. The deviations between the models 
Pi and Pi are only at the level of a few percent. This demonstrates that the entropy is 
indeed sensitive to the large-t behavior of p{t), while the short-time behavior of the models 
is irrelevant. 



13 



In the models P„, the spectral function decreases as t"". Since the functions fpnif) are 
related by derivatives with respect to t, their Fourier transform (28) can be calculated by a 
simple recursion relation given in Table I. Fig. 7 shows the spectral functions for n = 1 . . . 4 
along with the corresponding entropies. Plotting s as a function of 7 = 717 confirms again 
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FIG. 7: Results for the models Pn {n = 1...4) similar to Fig. 5. The left plot compares the spectral 
functions pn with a width parameter 7^ = 7/n and pL with 7 = 7 = 1; the right plot shows the 
entropies as a function of 7. 

the general expectation that the entropy is not determined by the small-t behavior, which 
is here Pn{t) ~ 1 — 7^, but rather by the large-time behavior of the spectral function. As 
expected, the faster p(t) decreases, the larger the entropy. Leaving the class of polynomial 
models, this trend is also obvious when considering the Gaussian model G in Fig. 8. 



IV. QCD 

In QCD, the entropy for various numbers of quark flavors, plotted as a function of T/Tc 
and scaled by the free limit, has a remarkably universal behavior as found in lattice calcula- 
tions [17]. I focus here on the representative case of the quenched limit of QCD, and point 
out briefly expected differences for the physical case. 

The gluon propagator has a transverse and a longitudinal part which leads to a cor- 
responding decomposition of the entropy. The contributions have the form (5) multiplied 
by the respective degeneracies [7]. The longitudinal modes are collective excitations whose 
spectral strength is exponentially suppressed for larger momenta, which leads to only a mi- 
nor contribution to the entropy. In the perturbative limit, it is of the order while the 
transverse modes yield a 0{g'^) term. Also for larger coupling the longitudinal contribution 
is rather small as demonstrated in the HTL calculations [7, 8], which actually might over- 
estimate the effect since it leads to a negative total entropy [35] at very large g. Taking 
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FIG. 8: Results for the Gaussian model G analogous to Fig. 5. 



therefore into account only the dominating transverse excitations with the propagator A, 
the resummed entropy reads 



2{Nl - 1) 1^ (^Imln(-A-^) + ImnReA^ . 



(30) 



In a self-consistent approximation, the propagator will have a residual gauge dependence, 
leading to an unphysical result for the entropy. However, as motivated before, a parameter- 
ization of the exact propagator by the dispersion relation and the width, which are gauge 
invariant, can be used in (30). 

The ansatz of the quasiparticle models [10] is to neglect the width altogether, and to 
describe the propagator simply by the perturbative asymptotic self-energy on the light-cone, 
which is a gauge invariant mass term. 



M2 = ^/t2, 
6 



(31) 



where = "i [36]. This 'minimal' resummation of the entropy can indeed nicely describe 
the lattice data for all temperatures above if an IR enhancement of the running coupling 
is permitted, for example in the form 



9\T) 



487r2 



IINM{\{T-Ts)/T,Y 



(32) 



For the physical number of degrees of freedom, dg = 2{N^ — 1), a fit of the parameters A 
and Tg/Tc leads in the entropy to small but systematic deviations from the lattice result for 
T > 2Tc, cf. Fit 1 in Fig. 9. This can be improved when considering dg as an additional 
fit parameter, which yields a value not too far from the physical one. While in the first 
reference [10] the longitudinal modes were considered as one possible explanation, it does 
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FIG. 9: The entropy of the SU(3) plasma in units of the free entropy. The symbols represent the 
lattice data [23]. Fits without width are depicted by the dotted and dashed line, for the latter dg 
was also fitted (see text). The full line is the Fit 3 with the width. The parameters are summarized 
in Table II. 





A 




dg 


c 


Fit 1 


10.5 


0.88 


16 (fixed) 




Fit 2 


5.2 


0.76 


17.5 




Fit 3 


2.6 


0.50 


16 (fixed) 


14.0 



TABLE II: The parameters of the fits shown in Fig. 9. 



not seem likely after the longitudinal contribution was found to be negative in the HTL 
calculation [7]. 

Taking now into account the width on the same footing as the mass (31), I will first 
consider a perturbative ansatz for 7. In the weak coupling limit, the width of a hard 
transverse gluon was obtained by Pisarski [22] in a resummed calculation, as 

0/1 \ '"-mag ^ ^'"'inag / / 

Several assumptions have been made here: (z) soft gluons are HTL dressed while intermediate 
hard gluons have a Lorentzian spectral function; (ii) the divergence from the static magnetic 
gluons is screened assuming that this sector of QCD can be parameterized by a mass mmag ~ 
g'^T. The first supposition follows the concept of a self-consistent calculation, hence 7 
appears also on the right hand side of Eq. (33) as a regulator next to the magnetic mass. 
Since the argument of the logarithm is basically (?~^, the width of moving excitations is 
enhanced compared to the width at rest, 7(0) ~ g^T. Implicit with (i) is the supposition 
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FIG. 10: The masses and the width according to the fits shown in Fig 9. 



7 = — 7^^ 111 7777^' (34) 



of a simple pole structure of the propagator, which is not warranted by any fundamental 
requirement. In fact, the result (33) could only be justified, to separate the pole from a 
branching point, if 7 <^ 'Timag (although physically the converse relation was considered 
more likely). While the constraint was necessary to explore details of the cut-off, the generic 
behavior 7 ~ g^\n.{l/g^)T is expected on general grounds [2]. To keep the connection to the 
result (33), I analyze the SU(3) entropy with the Lorentz spectral function and the width 
in the form [37] 

3 M2 , c 

rri 1 

where c parameterizes the soft cut-off. It is emphasized that because the functional relation 
between M and 7 is fixed, it is not obvious whether a fit is possible at all, or that introducing 
c in addition to the parameters in the coupling (32) will improve the result. The numerics 
shows, however, that this is indeed the case. The enhanced entropy, due to the width, nicely 
explains the small deviations of Fit 1 for T > 2 without having to adjust dg. As to be 
expected, the mass is somewhat larger than in the previous fits, cf. Fig. 10. More interesting 
is the behavior of the width near Tc. Because there s/sq is small, the mass and hence the 
coupling have to be large. At the same time, the width cannot be too large since it would 
over-compensate the decreasing effect of the mass. Within the ansatz (34) this implies that 
the logarithm has to become small. It is worth to emphasize that the optimal value of c, 
given in Table II, is surprisingly close to the value 

.-^^«13.7. (35) 

c 

so the width vanishes almost precisely at (and is indeed positive for all temperatures) 
[38]. Interestingly then, the condition 7 <^ w^mag) which was necessary in the derivation 
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of Eq. (33) (but originally considered to not represent the physical situation), can actually 
be fulfilled in a small vicinity, T+, of T^. Taking the next-to- leading logarithm result at 
face value leads to an estimate of the magnetic mass which is consistent below LIT^, where 
"^mag > 27. This estimate (including the 'predicted' range of applicability) is indeed in nice 
agreement with the lattice data [24] as shown in Fig. 11. The magnetic mass at Tc, which 



be 
cS 



o lattice 
^— T* estimate 
-- 0.27 

- 0.48^2 loop 
■• 27/ T 




FIG. 11: The lattice data [24] for the magnetic gluon mass, and the estimate based on (33) with 
the fitted width, which is meaningful only for T = T^, i. e., only for the left data point (see text). 
Also shown are the ^^T-fits with the 2-loop running coupling [24] (dash-dotted line), and with the 
adjusted coupling (32) from Fit 3 (dashed line). 

is difficult to calculate on the lattice, is estimated as 




(36) 



For larger temperatures, the behavior of the magnetic mass is consistent with the expectation 
"^mag = dg^T. Note that the fit from Ref. [24], which uses the 2-loop running coupling at 
the momentum scale 27rT, can be improved by using the coupling as obtained in Fit 3. 

It is emphasized that the magnetic and the electric screening masses have a qualitatively 
different large-coupling behavior: The electric screening mass, as known from lattice calcu- 
lations [24], becomes small near Tc. This is readily understood in the quasiparticle model 
[10], where 

dn 



TTIt 



noo(^ = 0,p-.0) = -2/iV, / — 

Jfc3 duj 



(37) 



Near T^,, the excitations are narrow (thus the quasiparticle picture is justified) and heavy, 
and the Debye mass 

-M/T j,2 (3g) 
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is Boltzmann-suppressed. This decrease cannot be expected, even tendentiously, from the 
next-to-leading order perturbative result [25], 



2 2 
"^D,n — ^D,0 



^/3iV; / 2mD,n 

1 -I m — 

27r \ mmag 



(39) 



Although apart from the more obvious solution, which is larger than uid q = (Nc/S^^'^gT, 
there is a second smaller solution (at Tg, with (36), it is mD,n ~ 0.4mD,o), the latter is 
unphysical because it does not approach the correct perturbative limit. 

Coming back to the discussion of the width, it is plausible from the properties of the 
entropy (e.g., from the fact SL{m — ^) — so mentioned in Sec. 3) that the width has to be 
rather small near Tc. However, the functional form (34), even as an extrapolation of the 
perturbative form similar to (31), is a priori not justified near Tc where M becomes large 
(see footnote [37]). Physically, one would here expect again a Boltzmann suppression of the 
heavy thermal fluctuations, hence 

^^Ae-^^g'^T. (40) 

For lack of better knowledge, 1 take this generic form as an ansatz for large coupling, and 
smoothly connect it to the 'perturbative' form (34) with the adjusted values of {A, Tg/T^ c}, 

by 

7* = (i-e)7 + e7, 

with 0(T) = arctan((r — T) /6) . Since the fit function s/so can basically be described 

by 3 parameters (say by the values at Tc and in the saturation-like regime, and by the slope 
at Tc), a conclusive determination of the parameters {A, 6, i/, T, 5} cannot be expected. In 
any case, within the enlarged parameter space the improvements in compared to Fit 3 
are only of the order of a few percent; the changes in the plot of the mass and the width, 
including the distinguished behavior at Tc, are almost invisible. This robustness of the results 
justifies a posteriori the usage of the perturbative ansatz (34) also for smaller temperatures. 
For larger T, after a distinguished maximum at 

T^^l.STc, (41) 

the ratio 7/T decreases very slowly. Different from what the parametric form of Eq. (34) 
might suggest, the width is even for rather large T to a good accuracy proportional to the 
mass; for T/Tc in [5,100], 

^ 0.69 - 0.02 In ^. (42) 

M Tc ^ ^ 

This underlines the fact that in this range of temperatures quasiparticle models can pro- 
vide only an effective description, while making a relation to the actual excitations may 
be difficult. Near Tc, on the other hand, the transverse hard excitations can be directly 
interpreted as quasiparticles as apparent from Fig. 12. This concept is beneficial since (up 
to rather large temperatures) the system is still strongly coupled: terms of higher order in 
g contribute significantly in the resummed entropy. 
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FIG. 12: The Lorentzian gluon spectral function from Fit 3 for T/T^. = 1.03, 1.35, 3 (u, k and 
p are in units of T). Shown here is the full phase space although the present approach can make 
statements only for hard momenta, of the order of T or larger. 

For definiteness, I have considered for QCD the case of a Lorentz spectral function. 
However, from the results of Sec. HIE it is obvious that the main result - a small width 
near Tc - should hold true also for other spectral functions, unless their Fourier transform 
has an exotic long-time behavior such as p{t) ~ which is not to be expected. 

V. CONCLUSIONS 

In this work it was shown for the deconfined SU(3) plasma that the collisional width 
(or damping rate) of hard gluons should be sizeable at intermediate temperatures, but 
has to become small near T^. While from an extrapolation of the parametric estimate 
7 ~ g'^\a{l/g)T (with the logarithm read as an enhancement factor) this result may seem 
surprising, a large width would be hard to reconcile with the small entropy near as 
established in lattice calculations. For QCD with dynamical quarks, the rescaled entropy 
has a similar temperature dependence [17]. Therefore, the principal result of a small width 
near the transition carries over from the quenched to the physical case. Although in full 
QCD s/sq is slightly larger at Tc, which is related to the nature of the transition, and the 
widths of the hard transverse gluons and quark particle-excitations might not fully vanish, 
a quasiparticle picture of the strongly coupled QCD plasma close to Tc appears justified. 

There are several interesting implications of the characteristic temperature dependence of 
the damping rate besides those for the screening properties discussed above. As the inverse 
of the mean free path A, it is closely related to transport properties as, e.g., equilibration 
times. Expecting a critical slowing down, this link provides another indication that the 
width has to become small near T^.. Another quantity, which is of particular interest with 
regard to the interpretation of SPS and RHIC experiments, is the radiative energy loss of 
hard quarks and gluons transversing the plasma. The results derived in [26] for a system of 
length L under the assumption of (several) independent scatterings, i. e. mo ^ 7, which is 
not unrealistic in the situation of interest, see below, are characterized by the energy scale 



(43) 
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In the Landau- Pomeranchuk-Migdal regime, for parton energies E > E^r, the energy loss 
reads [26] 

-AE = ICr a-fmlL'^ In 7L , (44) 

where R indicates the color representation of the parton. As argued in [27], the expected 
'critical' behavior of the screening mass would lead to a reduced energy loss at tempera- 
tures near Tc, possibly explaining the absence of jet quenching at SPS energies [28]. With 
the temperature dependence of the width suggested here, this effect would be even more 
pronounced. At the same time, the adjustment of the parameters to lattice data may allow 
for a realistic estimate. In order to describe the Debye mass in a simple way, without any 
assumptions, I make use of the empirical observation that for the relevant temperatures 

niD ^ 2.77 , (45) 

cf. Fig. 13. The resulting behavior of the energy scale E^r, shown in Fig. 14, changes 




FIG. 13: The lattice data [24] for the gluon Debye mass and the rescaled width from Fit 3. 



drastically at 

Te^I.ST,. (46) 

Close to Tc, the energy E^r becomes very small. Already for slightly larger temperatures it is 
well above the scale -Elpm = '^d/t which is relevant for the energy loss in the Bethe-Heitler 
regime. A similar sudden onset, also at Te, is found for the medium induced energy loss, cf. 
Fig. 15. For larger temperatures the estimate agrees basically with other results, while close 
to Tc it becomes very small and would be hard to observe experimentally. Similar small 
results have been obtained for the corresponding SPS energies in Ref. [29]. 

In summary, it has been argued from the reduced number of degrees of freedom near 
the transition temperature that the width of hard excitations has to become small near 
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FIG. 14: The temperature dependence of the energy scales E^r and -Blpm for the radiative energy- 
loss. The length of the medium was set to L = 5fm, and Tc = 170 MeV. 

Tc- While for QCD this was demonstrated under the assumption of a Lorentz spectral 
function, the propagator may have a more complicated pole structure. For hot QED (where 
no magnetic mass exists), the fermion propagator has been calculated in an Bloch-Nordsieck 
approach [30]. The result is infrared finite, and an entire function of the energy. The spectral 
function is nonetheless strongly peaked, with a characteristic width ~ e^ln(l/e)T. Since 
in Fourier space it decreases faster than an exponential, the effect for the entropy will be 
even more pronounced than for a Lorentzian. While the situation may be different in QCD, 
a spectral function with p{t) ~ which has little effect on the entropy, seems hard to 
explain. Although from thermal properties obtained in lattice QCD little can be inferred 
about the analytic structure of the propagator, the general result of small widths near is 
arguably robust. This implies a characteristic change of several phenomenologically relevant 
quantities at T* ^ ^ Te ^ LST^. 
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APPENDIX A: SOME PROPERTIES OF AS 

In the following it is argued that the entropy is generally increased for a non-zero width, 
i. e. As > 0. Furthermore, considering the Lorentz spectral function, the expansion of 
•Sl('^5 7) is calculated for small arguments. Finally it is proved that S2,(m, 7) is for m = 7 
equal to the Stefan-Boltzmann value Sq = s^'^\m = 0) = ^T^. 
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FIG. 15: The radiative energy loss from Eq. (44) for L = 5fm, unless indicated otherwise. The 
symbol is an estimate from Ref. [26]. The dashed line results from a 2-loop running coupling, 
with mD set to mD,o and a constant 7 to match the result from [26]. The full line is the estimate 
with the T-dependent width and the adjusted coupling, as is the dotted line where L = 3.5 fm was 
assumed to illustrate the scaling behavior. 

In order to prove (under the assumption of a unique dispersion relation ujk) that As > 0, 
consider the relevant integral in (10) in the form 



where / = arctanA — A/(l + A^), and A = ImA/ReA. Changing the integration variable to 



the integrand becomes a product of three factors, of which /(A) is an odd function. The sign 
of / is, thus, determined only by the other two terms, which can be discussed on the basis of 
the relation co'(A), whose inverse is shown in Fig. 2. In particular, oj — > {0, uJk — 0, 0;^ + 0, 00} 
corresponds to A — > {O"*", +cxo, — cxd, 0~}. Then, via the inverse derivative dX/doj and the 
properties of the propagator listed in Section 2, it is easily inferred that du/dX is positive 
and that it vanishes for |A| —>■ 00, is non-zero for A ^ 0^, and diverges for A ^ 0^, cf. 
Fig. 16. The product du/dXdn/dT, however, vanishes for small negative A due to the 
second factor, which for the corresponding large u is exponentially suppressed. At small 
positive A, on the other hand, this factor is Bose enhanced, dn/dT ~ cu"^, while it goes to 
the value dn/dT\^^ for A — > ±00. Therefore, the integral of the product {du/dX) (dn/dT) 
over [—00, 0] is finite, while for the interval [e, +00] it approaches -|-cxo for e ^ 0. This shows 
that the integral J, whose integrand includes the odd function /(A), is positive (and finite 
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FIG. 16: Two functions used in the argumentation that As > 0, for the case of the propagators 
from Fig. 1. 



due to / ~ for small A), hence As > 0. 

Turning now to the expansion of the entropy ^^(m, 7) = s^^\m)+AsL{m, 7) for a Lorentzian 
spectral function, I introduce the notation T = 27 and follow the remark below Eq. (10), 
writing 

where h = arctan(ru;/(u;^ — cu^)). Considering first the case m = 0, i.e., Um = k, the two 
terms of the fc-integrand decrease as k~'^. In the subtracted integral 



1= I i^-^ 



UJ 



fc3 \dT 27r-Jo 



dk k' 



k^-co^ 



the substitutions k = xuo and a = T /uo lead to 

,2 ^no / x2 — 1 



J 



dx X 



[x^ — 1)2 + a2 X 



27r2 7o 



a;2 TT / xl/2 ujV 

— j= (Vl + a2-l) = 

27r2 2^/2 ^ ^ Stt 



The remaining a;-integral yields 

f°° duo dn ^ T ^ 

J = / — -7^1 = — r + ... 

Jo 71 dT 24 

From AsL^m, T) = jdV J — VJ, and since Asiim, F = 0) = (for any m), one obtains 

Asz.(m = 0,F) = ^F2 + ... 
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Furthermore, it is obvious that derivatives of any order of Asx,(m, F) with respect to m 
vanish at F = 0. Therefore, the leading term in m in the expansion of the total entropy 
comes entirely from the contribution s^^\m), which is well known, and one arrives at the 
expression (19). 

Finally, the fact that SL{m,'-f = m) = Sq holds not only in the limit of small m and 7 
is readily proven by verifying dsL{m,'~f = m)/dm = 0. After taking the derivative of the 
integrand of the total entropy [39], the /c-integration yields zero, indeed. As an aside it is 
mentioned that for 7 > m, the poles of the propagator (16) are purely imaginary for some 
range of momenta, see footnote [34] . 
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